home *** CD-ROM | disk | FTP | other *** search
/ MacFormat 1995 April / macformat-023.iso / Shareware City / Developers / RLaB / rlib / std.r < prev    next >
Encoding:
Text File  |  1994-09-23  |  733 b   |  29 lines  |  [TEXT/RLAB]
  1. //-------------------------------------------------------------------//
  2.  
  3. //  Syntax:    std ( A )
  4.  
  5. //  Description:
  6.  
  7. //  Compute the Standard Deviation. For a vector (row-matrix), std()
  8. //  returns the standard deviation. For a MxN matrix, std() returns a
  9. //  row-matrix containing the standard deviation of each column.
  10.  
  11. //  See Also: mean, rand, srand
  12. //-------------------------------------------------------------------//
  13.  
  14. std = function(x)
  15. {
  16.   if(class(x) != "num") { error("std() requires NUMERICAL input"); }
  17.  
  18.   m = x.nr;
  19.   if( m == 1 ) 
  20.   { 
  21.     return sqrt( sum( (x - mean(x)) .^ 2 ) / (x.nc - 1) );
  22.   else
  23.     for( i in 1:x.nc) {
  24.       s[i] = sqrt( sum( (x[;i] - mean(x[;i])) .^ 2 ) / (x.nr - 1) );
  25.     }
  26.     return s;
  27.   }
  28. };
  29.